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FOREWORD 


This  report  covers  work  carried  out  at  AFAPL'  s Turbo  Structures 
Research  Laboratory  (TSRL)  on  the  transient  structural  response  of  an 
isotropic  cantilevered  plate  subjected  to  normal  impact  by  a ballistic 
pendulum.  The  effort  was  intended  as  a vehicle  for  evaluating  the 
methods  of  pulsed  laser  holography  and  finite  element  analysis  as  they 
relate  to  the  study  of  transient  structural  dynamics.  This  is  an  area 
which  bears  directly  on  the  problem  of  foreign  object  damage  to  turbine 
engine  components. 

The  program  was  a combined  expermental/analytical  effort.  The 
experimental  portion  utilized  a pulsed  ruby  laser  to  obtain  holographic 
interferograms  of  the  plate's  deformation  following  impact.  The 
analytical  portion  of  the  work  consisted  of  mathematically  modelling 
the  plate  using  finite  element  techniques  and  studying  the  model's 
response  to  impact  using  the  general  purpose  finite  element  program, 

NAS TRAN. 

The  work  was  performed  in  the  Turbine  Engine  Division  of  the  Air 
Force  Aero  Propulsion  Laboratory,  Air  Force  Systems  Command,  Wright- 
Patterson  Air  Force  Base,  Ohio  under  Project  3066,  Task  12,  and  Work 
Unit  21.  The  effort  was  conducted  by  Dr.  James  C.  MacBain  of  the 
Propulsion  Branch. 

The  author  is  indebted  to  Mr.  Bruce  Tavner  for  his  very  competent 
technical  assistance  in  the  laboratory  and  to  Miss  Helen  Davis  for 
typing  the  manuscript. 
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SECTION  I 
INTRODUCTION 


I 


I 


Foreign  object  damage  in  both  military  and  civilian  turbine 
engines  is  a seldom  but  serious  problem  both  in  terms  of  cost  and 
safety.  For  example,  the  Air  Force  Inspection  and  Safety  Center 
states  that  2816  birdstrikes  have  been  reported  Air  Force-wide  be- 
tween 1966  and  1973.  These  birdstrikes  resulted  in  a loss  of  7 lives, 
14  aircraft,  and  a cost  to  the  Air  Force  of  $74  million  dollars. 

These  figures  serve  to  underline  the  fact  that  there  is  a definite 
need  for  both  basic  and  applied  research  in  the  area  of  impact- 
tolerant  turbine  engine  blading.  More  sophisticated  design  tools 
and  impact  theories  are  required.  In  this  vein,  this  report  covers 
work  carried  out  at  AFAPL's  Turbo  Structures  Research  Laboratory 
(TSRL)  on  the  transient  structural  response  of  an  isotropic  canti- 
levered plate  subjected  to  normal  impact  by  a ballistic  pendulum. 

The  program  was  a combined  experimental/analytical  effort.  The 
experimental  portion  utilized  a pulsed  ruby  laser  to  obtain  holo- 
graphic interf erograms  of  the  plate's  deformation  following  impact. 

The  analytical  portion  of  the  work  consisted  of  mathematically 
modelling  the  plate  using  finite  element  techniques  and  studying  the 
model's  response  to  impact  using  the  general  purpose  finite  element 
program,  NASTRAN. 

The  specific  aims  of  this  research  effort  were  threefold.  First, 
the  advantages  and  disadvantages  of  using  pulsed  laser  holography  for 
transient  structural  analysis  were  studied.  Second,  it  was  intended 
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to  use  the  experimental  pulsed  holography  results  as  a verification 
of  results  obtained  analytically  using  the  computer  program,  NASTRAN. 
Finally,  the  research  effort  served  to  provide  knowledge  of  and 
experience  with  pulsed  laser  holography  - information  that  will  be 


ncpfnl  in  fnfiirp  TQRT  rocnordi  of fnrfn 
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SUMMARY  OF  RESULTS 


An  aluminum  cantilevered  plate  measuring  3"x7"x.l875"  was  struck 


with  a ballistic  pendulum  consisting  of  a .65"  diameter  steel  ball 


attached  to  a wire  on  a pivot.  The  resulting  plate  response  was 


measured  for  a specified  time  after  impact  by  double  exposure  holographic 


interferometry  using  a pulsed  ruby  laser.  The  plate's  normal  displace 


ment  was  experimentally  determined  for  times  after  impact  ranging  from 


2 to  33ps.  Photographs  of  the  doul !e  exposure  holograms  from  these 


tests  are  shown  in  Figures  4. 1-4. 5.  The  normal  displacement  based  on 


four  of  the  test  runs  (times  after  impact  of  4,  6,  12,  and  18ps)  is 


shown  as  a function  of  plate  geometry  in  Figures  4.10-4.13 


The  flexural  wave  velocity  was  computed  from  the  holographic 


interf erograms  by  plotting  the  plate  wave  position  versus  time  after 


impact  and  was  found  to  be  Cf  = .102  in/ps.  This  is  in  good  agreement 


with  the  theoretical  Rayleigh  surface  wave  velocity  of  .112  in/ps 


a difference  of  8.5% 


A parallel  numerical  study  was  conducted  using  the  finite  element 


computer  program  NASTRAN  to  compute  the  cantilevered  plate's  transient 


response.  The  results  were  in  good  agreement  with  the  experimental 


findings.  The  plate's  normal  displacement  based  on  finite  element 


analysis  is  shown  plotted  in  Figures  4.10-4.13  as  a function  of  plate 


geometry  (dashed  lines).  Contour  plots  of  the  plate's  normal  displace 


ment  for  different  times  after  impact  were  also  generated  by  NASTRAN 


and  are  shown  in  Figures  5. 3-5.5 
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The  study  demonstrated  the  feasibility  and  utility  of  using 
pulsed  laser  holography  to  study  transient  structural  response. 
In  addition,  it  provided  increased  confidence  and  experience  in 
the  use  of  NASTRAN's  transient  analysis  capability. 
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SECTION  III 


EXPERIMENTAL  SET  UP  AND  PROCEDURE 


3.1  Physical  Configuration 

The  test  piece  for  the  experimental  portion  of  the  impact  analysis 
program  was  a 6061-T6  aluminum  plate  measuring  12"  in  length,  3"  in 
width,  and  3/16"  thick.  The  plate  was  fixed  between  two  steel  blocks 
having  a total  weight  of  33  lbf  such  that  it  was  cantilevered  and  had 
a free  length  of  7".  The  weight  of  the  cantilevered  portion  of  the 
plate  was  .394  lbf.  The  plate  and  jig  are  shown  in  Figure  3.1. 

The  plate  was  impacted  normally  by  a steel  ball  weighing  .043  lbf 
at  a point  lying  on  its  long  axis  and  located  3"  above  its  fixed  end  as 
shown  in  Figure  3.1.  The  ball  was  soldered  to  a thin  wire  that  in  turn 
was  fixed  to  a pivot  located  a distance  above  the  impact  point  forming 
what  is  known  as  a ballistic  pendulum.  The  impact  sequence  was  initiated 
by  suspending  the  steel  ball  from  an  electromagnet  which  was  then  switched 
off  allowing  the  ball  to  swing  down  and  strike  the  plate.  Just  prior  to 
impact,  the  ball  interrupted  a continuous  wave  laser  beam  passing  behind 
the  plate  (see  Figure  3.2)  causing  a photo  diode  to  transmit  a 10V 
signal  that  initiated  the  pulsed  laser  firing  sequence.  The  timing  and 
electronics  involved  in  the  pulsed  laser  firing  sequence  will  be 
addressed  in  more  detail  in  a later  section. 


The  optical  set  up  for  making  the  hologram  is  also  shown  in  Figure  3.2. 

The  placement  of  the  optics  is  typical  of  that  used  to  make  transmission 
holograms  with  an  off-axis  holographic  set  up,  and  the  reader  is  referred 

I 

5 


mam m 


m Bflu 


. a # ,-•  ... 


Figure  3.1  - Cantilever  Plate  and  Pendulum  Geometry 


Experimental  Set-up  for  Making  Holographic  Interferograms 


to  Reference  1 for  details.  The  holograms  were  made  using  a pulsed 
ruby  laser  (Apollo  Model  22HD) . The  laser  puts  out  2.5  joules/pulse, 
20-50ns  In  width,  which  necessitates  the  use  of  dielectric  mirrors 
and  doubly  concaved  lenses  (F.L.  = -40mm)  when  working  with  the  un- 
expanded laser  beam.  While  TSRL's  pulse  ruby  laser  is  capable  of 
rapid  double  pulsing  down  to  pulse  separations  of  lps,  this  was  not 
done  in  the  present  experiment.  Because  of  timing  constraints  placed 
on  the  ruby  laser  (see  Section  3.2),  the  holographic  interf erograms 
were  made  by  first  exposing  the  photographic  plate  to  a single  manually 
initiated  laser  pulse  prior  to  impact  and  then  again  exposing  the 
photographic  plate  some  time.  At,  after  impact.  This  second  exposure 
was  also  done  in  the  single  pulse  mode  but  in  this  case  was  triggered 
automatically.  This  will  be  discussed  in  greater  detail  in  the  next 
section. 

3.2  Electronic  Timing  Circuitry  for  Firing  Laser 

When  the  ballistic  pendulum  impacts  the  cantilever  plate,  flexural 

waves  spread  out  from  the  impact  point  at  a velocity  approaching  the 

shear  wave  velocity,  Cg  = £,  of  the  material  (Ref.  2).  For  6061-T6 

P 

aluminum,  the  flexural  waves  will  have  a phase  velocity  approaching 
.114  in/us.  This  means  that  when  the  plate  is  'impacted,  the  resulting 
flexural  waves  will  only  take  about  16us  to  travel  the  1.5  inches  to 
the  free  sides  of  plate.  Now,  one  of  the  objectives  of  the  present 
study  was  to  analyze  the  plate  deformation  just  after  impact,  i.e., 
at  times  after  impact  prior  to  significant  wave  reflection  off  the 


plate's  boundaries.  For  acceptable  results,  this  places  an  upper 
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bound  on  the  time  for  firing  the  second  laser  pulse  of  about  30ys 


(wave  reflection  will  occur  off  the  closest  free  edges  of  the 


plate  at  t = 16ps).  This  in  turn  places  some  constraints  on  the 


firing  sequence  for  the  ruby  laser  since,  as  shown  in  Figure  3.3,  the 


ruby  laser  can  fire  only  after  200ps  have  elapsed  - the  amount  of 


time  necessary  for  the  flash  lamps  to  energize.  This  is  true  for 


mo  ie  of  operation.  Hence,  in  order  for  the  laser  to  lase  automatically 


at  some  time  after  impact  in  the  0-30ps  range,  it  would  have  to  be 


triggered  (signal  sent  to  Master  Sync)  at  some  time  prior  to  impact 


This  was  done  using  a timing  delay  oscilloscope  (Tektronix  535)  as 


follows 


Prior  to  making  the  double  exposure  hologram,  the  time  that  it 


took  the  ball  to  go  from  the  trigger  laser  beam  to  the  plate  was 


measured  by  an  electronic  counter.  This  is  shown  as  T._  in  Figure  3.4 


This  distance  traversed  by  the  ball  in  going  from  A to  C was  approximately 


1 inch  and  T.„  was  typically  in  the  .02  second  range  (20,000ps).  Upon 


specifying  the  time  after  plate  impact  to  be  observed,  T 


the  fact  that  Pulse  #2  of  the  laser  fires  at  lOOOps  into  the  firing 


sequence,  the  delay  time,  T 


was  determined  from  the  relationship 


lOOOps 


The  delay  oscilloscope  was  then  set  so  that  a signal  was  sent  to  trigger 


the  ruby  laser  after  a delay  of  T.  seconds.  The  Pockels  cell  voltage 


for  Pulse  //I  of  the  ruby  laser  was  set  to  zero,  thus  eliminating  Pulse  #1 


Figure  3.3  - Timing  Sequence  for  Pulsed  Ruby  Laser 
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With  the  steps  mentioned  above  carried  out,  the  actual  impact 
test  was  run.  First,  the  magnet  was  energized  and  the  pendulum 
attached  to  it.  The  first  exposure  of  the  stationary  cantilever  plate 
was  then  made  by  manually  firing  the  ruby  laser  (Pulse  if 2,  only). 

The  ruby  laser  then  was  energized  again,  the  magnet  was  switched  off, 
allowing  the  ballistic  pendulum  to  swing  down,  interrupt  the  trigger 
laser  beam  (point  A),  initiate  the  ruby  laser  firing  sequence  (point  B) , 
and  impact  the  plate  (point  C).  At  T^pj,  seconds  after  initial  contact 
between  the  ball  and  plate.  Pulse  if 2 of  the  ruby  laser  fired,  exposing 
the  photographic  plate  for  the  second  time.  During  this  impact  sequence, 
the  time  T was  again  measured  as  it  varied  somewhat  between  tests. 
Knowing  the  value  of  T for  the  actual  test,  the  actual  value  of  the 
time  after  impact,  T^^p^  could  be  computed  from  equation  (3.1).  A 
schematic  of  the  electronics  used  in  the  experiment  is  shown  in 
Figure  3.5.  A switch  was  put  in  the  electronic  counter  portion  of  the 
circuit  (located  at  bottom  of  Figure  3.5)  in  order  to  facilitate  quick 
monitoring  of  the  present  scope  delay  time,  T , or  the  time  it  took 
the  ball  to  go  from  the  trigger  laser  to  the  plate,  T^c 

3.3  Hologram  Processing 

The  photographic  plates  used  in  the  experiment  were  4"x5"  Agfa- 
Gevaert  10E75.  The  photo  cell  and  storage  oscilloscope  shown  behind 
the  hologram  in  Figure  3.5  were  used  to  measure  the  reference  to  object 
beam  ratio  and  pulse  amplitude  of  the  holographic  set  up.  The  beam  ratio 
for  the  tests  was  approximately  2 to  1 (reference  beam  to  object  beam). 
The  combination  of  photo  cell  and  storage  oscilloscope  could  only  give 


Figure  3.5  - Electronic  Circuitry  for  Making  a Double  Exposure  Hologram 
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relative  light  intensity  measurements,  i.e.,  the  amplitude  of  the 
reference  beam  relative  to  the  object  beam  or  the  relative  amplitude 
between  the  first  and  second  plate  exposures.  This  sufficed,  however, 
for  the  present  experiment.  With  the  photo  cell  placed  1"  behind  the 
hologram  plate  holder,  "acceptable"  holograms  resulted  for  oscilloscope 
readings  of  .4  volts  and  .6  volts  with  and  without  a photographic  plate 
in  the  holder,  respectively.  The  resulting  holograms  tended  to  be  dark 
but  this  was  corrected  by  bleaching  the  developed  photographic  plates. 
The  developing  process  was  carried  out  as  follows: 

1.  4 minutes  in  Kodak  D-19  developer; 

2.  30  seconds  in  stop  bath; 

3.  3 minutes  in  fixer; 

4.  Dry  with  blow  dryer; 

5.  Bleach  in  Potasium  Ferricynanide  (15g  of 
K^Fe(CN)^  in  1000  ml  distilled  water)  - agitate 

until  plate  becomes  milky  white,  approximately 
1 minute. 
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SECTION  IV 


EXPERIMENTAL  DATA  REDUCTION  AND  RESULTS 


Figures  4.1  through  4.5  show  the  cantilevered  plate's 
response  for  times  after  impact  ranging  from  2ys  to  33ys.  For  the 
sake  of  increased  fringe  clarity,  the  photographs  taken  of  the 
holograms  show  the  immediate  area  of  impact.  Recall  that  the 
impact  point  lies  on  the  vertical  axis  of  the  plate  and  3 inches 
above  the  clamped  base.  The  numbers  in  parentheses  refer  to  the 
raw  data  numbering  system  used  to  denote  each  test  run.  The 
6061-T6  aluminum  plate  was  painted  with  a flat  white  paint  and 
a grid  having  1"  by  .5"  increments  was  scribed  on  it.  In  addition, 
1/8"  increments  were  scribed  on  the  plate  centerline  from  the 
impact  point  to  1"  below  it. 

Referring  to  Figures  4.1  through  4.5,  the  fringes  concentrically 
located  around  the  impact  point  (ball  impacting  from  rear)  represent 
loci  of  constant  displacement  on  the  plate.  The  fringes  are  a contour 
map  of  the  flexural  waves  caused  by  the  impact.  They  travel  outward 
with  time  until  they  reach  the  free  edge  of  the  plate  at  about  T=12ps 
and  are  reflected.  It  is  evident  from  the  photographs  that  either 
the  timing  measurements  are  in  error  by  as  much  as  ± 2ps  or  that  there 
was  some  scatter  in  the  magnitude  and  duration  of  the  impulse  load 
imparted  to  the  plate  by  the  steel  ball  pendulum.  The  latter  reason 
is  thought  to  be  the  case  because  of  a permanent  magnetic  field  that 
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Figure  4.5  - Double  Exposure  Hologram  of  Cantilever  Plate  (28,  30,  and  33ps) 
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was  induced  in  the  electromagnet  used  to  release  the  steel  ball. 

The  permanent  magnetic  field  varied  somewhat  throughout  the  dura- 
tion of  the  impact  tests.  The  magnetic  force  field  would  have  the 
effect  of  slowing  up  the  ball's  release  and,  consequently,  decreasing 
its  impact  force.  This  fact  is  corroborated  by  Figures  4.6  and  4.7 
where  the  force  time  profile  is  seen  for  four  successive  ball  impacts 
on  a pressure  transducer  (Piezotronics , PCB  118A-B,  2ps  rise  time) 


1 


temporarily  imbedded  in  the  plate  at  the  impact  point.  The  pressure 
amplitude  (proportional  to  force)  in  test  1 is  20%  less  than  the  tests 
2 through  4.  While  future  work  in  this  area  should  strive  for  a more 
repeatable  loading  system,  the  present  data  scatter  was  found  to  be 
acceptable  within  the  scope  of  the  present  research  effort.  The  load 
profile  shown  in  Figures  4.6  and  4.7  was  also  utilized  to  determine 
the  loading  input  for  the  finite  element  analysis  (Section  5.1). 

The  fringe  pattern  position  and  the  corresponding  times  after 
impact  in  Figures  4. 1-4. 5 can  be  used  to  compute  the  velocity  of  the 
flexural  waves  in  the  plate.  Using  the  outermost  visible  fringe  for 
times  after  impact  ranging  from  2ys  to  13ps,  one  can  get  a plot  of 
wave  position  as  a function  of  time.  The  slope  of  this  curve  will  be 
the  flexural  wave  velocity.  Using  the  fringe  pattern  positions  and 
corresponding  times  in  Figures  4.1  through  4.3,  the  curve  shown  in 
Figure  4.8  was  generated.  Because  of  the  aforementioned  scatter  due 
to  load  repeatability,  a linear  least  squares  routine  was  used  to 
generate  the  curve  through  the  data  points.  The  slope  of  the  curve 
gives  the  flexural  wave  velocity  as  C^.  = .1024  in/us. 
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tion  of  Time  After  Impact 
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It  is  of  interest  to  compare  the  experimental  value  of  the 


flexural  wave  velocity  to  that  obtained  using  the  three-dimensional 


equations  of  elasticity  (Ref.  2,  3,  and  A).  According  to  classical 


plate  theory,  the  flexural  wave  velocity  in  an  infinite  plate  approaches 


that  of  a Rayleigh  surface  wave  when  the  wavelength.  A,  becomes  small 


compared  with  the  plate  thickness,  h.  In  the  case  of  the  present 


impact  tests,  the  steel  ball  made  contact  with  the  plate  over  a 


circular  area  of  less  than  .03125  inches  in  diameter.  Now,  the  wave 


length,  A , is  a function  of  the  pulse  shape  and  plate  contact  area 


and  if  the  pulse  loading  is  broken  into  its  Fourier  components,  the 


largest  component's  wavelength  will  be  on  the  order  of  twice  the 


diameter  of  the  contact  area.  Hence 


the  plate  with  h = .1875  is  h/A  = 3.  All  other  values  of  h/A  for  the 


higher  Fourier  components  will  be  greater  than  this  value.  For  values 


of  h/A  > 3,  the  flexural  wave  velocity  in  a plate  closely  approximates 


that  of  a Rayleigh  surface  wave  and  can  be  obtained  from  the  expression 


is  the  shear  wave  velocity 


where  C 


lbm.  Placing  these  values  into  equation  (A.l) 


This  is  within  + 8.5%  of  the 


experimental  value 
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The  magnitude  of  the  plate's  displacement  can  be  obtained  from 
the  fringe  photographs  using  the  expression  (Ref.  5 and  bj: 

6.(n  + „ ) = for  N = ±1,  ±2,  ±3  ....  (4.2) 

o v 2 

where  j - displacement  vector 

- unit  vector  in  direction  from  object 
to  illumination  source  (object  beam) 

n - unit  vector  in  direction  from  viewer 
v 

(through  hologram)  to  object 

A - wavelength  of  laser  used  to  make  the 
hologram 

N - fringe  order 

The  vectors  given  by  equation  (4.2)  are  shown  in  the  context  of  the 
present  experimental  geometry  in  Figure  4.9. 


Carrying  out  the  dot  product  in  equation  (4.2)  yields: 


6 1 n + n Cos (n  , n ) = 
1 1 o v1  o v 


(2N+  1)  A 
2 


(4.3) 


letting  | 6 | =6  and  using  Figure  4.9  gives 


n + n = 

1 n v 1 


1.975  i + .222  j | = 1.987 


(4.4 
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Hence,  the  displacement,  6,  in  the  direction  that  bisects 

n and  n is: 
o v 


6 = 


(2N + 1)A 

(2) ( .973) (1 . 987) 


(2N  + 1)A 
3.874 


(4.6) 


The  displacement  normal  (perpendicular)  to  the  plate  surface  differs 
from  that  given  in  equation  (4.5)  by  only  Cos  f^Q,  n \=  -994. 


h I 

i p 


Dividing  equation  (4.6)  by  .994  gives 

5 _ (2N+  1) A , for  N = ±1,  ±2,  ±3  ...  (4.7) 

n 3.853 

Equation  (4.7)  will  be  used  below  to  compute  the  values  of  the 
plate's  normal  displacement. 

Figures  4.10-4.13  show  the  plate  displacement  as  a function  of 
the  distance  from  the  impact  point  along  a line  from  the  impact  point 
to  the  free  edge  of  the  plate  (Section  A-A  in  the  figures)  for  times 
after  impact  of  4,  6,  12,  and  18ps.  Also  shown  in  these  figures  are 
the  normal  displacement  curves  based  on  finite  element  analysis  which 
will  be  discussed  in  Section  V, 

Polaroid  type  55  P/N  film  was  used  to  photograph  the  double 
exposure  holograms  using  a 4"x5"  format  view  camera  with  a Polaroid 
film  holder.  Polaroid  55  P/N  film  produces  both  a positive  and  negative 
print.  The  negative  print  was  placed  in  a standard  photographic 
enlarger  such  that  the  enlarged  view  of  the  displacement  fringes 
could  be  used  to  more  accurately  determine  the  corresponding  dis- 
placement. Using  the  enlarger,  fringes  as  high  as  N=70  could  be 

28 


V 

r! 

■4. 


(Nirf) 


AFAPL-TR-7b-56 


Figure  4.10  - Normal  Displacement,  6 . Vs  Distance  From 

n 

Impact  Point,  Y,  at  T=4ps 
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observed.  (Note  that  this  technique  could  be  conveniently  expanded 
to  include  a digital  X-Y  plotting  table  as  the  projecting  surface 
for  the  enlarger.) 

In  contrast  to  time-average  holography  where  the  white  fringe 
denoting  zero  displacement  is  the  most  intense  and  stands  out  clearly, 
the  fringe  representing  zero  displacement  in  double  exposure  holography 
is  of  the  same  intensity  as  its  neighboring  higher  order  fringes. 

Hence,  in  order  to  get  a quantitative  plot  of  the  displacement  based 
on  a single  hologram,  some  a priori  information  on  the  plate's 
displacement  response  is  necessary.  Referring  to  Figures  4.1,  4.2, 
and  4.3,  the  maximum  plate  displacement  occurs  at  the  impact  point 
and  decreases  in  magnitude  as  one  travels  outward  from  it.  Since 
the  flexural  displacement  is  a wave,  the  positive  displacement  at 
the  impact  point  will  be  followed  by  a smaller  negative  one  at  some 
distance  from  the  impact  point  and  then,  as  one  travels  further  away 
from  the  impact  center,  a return  to  the  undeformed,  zero  displacement, 
portion  of  the  plate.  The  peak,  negative  displacement  manifests 
itself  by  a widening  of  the  fringes  where  the  slope  of  the  negative 
displacement  goes  to  zero.  This  can  be  seen  in  Figure  4.3  in  the 
12ps  (107)  photograph.  Using  this  type  of  reasoning  the  plate's 
normal  displacement  could  be  plotted.  While  this  approach  was  found 
to  be  sufficient  for  the  simple  deformation  pattern  that  the  plate 
experienced,  more  complex  displacements  would  necessitate  more 
sophisticated  approaches  such  as  multiple  double  exposure  holograms. 

As  an  aside,  another  approach  that  could  be  used  to  determine  the 
zero  displacement  fringe  would  be  to  use  an  optical  bench  telescope 
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to  scan  across  the  image  of  the  hologram.  By  scanning  along  a 
horizontal  line  passing  through  the  impact  point,  the  concen- 
trically located  fringes  will  appear  to  converge  toward  one  of 
the  circular  fringes.  In  other  words,  the  impact  point  will 
appear  to  be  a fringe  "source"  with  fringes  traveling  toward 
the  one  stationary  circular  fringe.  Fringes  located  outside 
the  stationary  circular  fringe  will  appear  to  travel  in  the 
direction  of  the  impact  point  and  the  stationary  circular  fringe. 
The  stationary  fringe  is  a loci  of  zero  displacement  and  by 
viewing  the  convergence  (traveling)  characteristics  of  the 
adjacent  fringes,  it  can  be  located. 

Figure  4.14  shows  a plot  of  the  plate's  normal  displacement 
along  its  free  edge  (Section  B-B)  at  T=24ps  after  impact.  The 
displacement  curve  is  based  on  the  third  fringe  photograph  shown 


in  Figure  4.14  and  shown  enlarged  in  Figure  4.15.  The  fringes 
are  numbered  in  Figure  4.15.  Along  the  free  edge,  the  undeformed 

I plate  was  used  as  the  zero  reference  point.  In  addition,  fringes 

that  curved  from  one  point  on  the  free  edge  to  another  provided  a 


; < convenient  indicator  of  points  of  equal  displacement  on  the 

$ 
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X (in),  SECTION  B-B 


Figure  4.14  - Normal  Displacement  Along  Plate  Free  Edge 
at  T=24ps 
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Figure  4 . 15 


Enlarged  View  of  Impacted  Plate  at  T=24ys  Showing 
Numbered  Fringes 
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SECTION  V 

FINITE  ELEMENT  ANALYSIS  AND  RESULTS 
5.1  Finite  Element  Model 

The  numerical  portion  of  the  study  was  based  on  a finite 
element  analysis  of  the  cantilever  plate  using  the  general  purpose 
finite  element  computer  program,  NASTRAN  (Navy  Nastran,  Level 
15.2.0).  The  model  geometry  and  orientation  for  the  cantilever 
plate  is  shown  in  Figure  5.1.  The  mesh  consists  of  304  nodes  that 
connect  the  165  quadrilateral  plate  elements  (CQUAD2) . The  nodes 
at  the  base  of  the  plate  were  fixed  against  both  translation  and 
rotation  to  simulate  the  cantilever  condition.  Because  the  impact 
load  acts  symmetrically  with  respect  to  the  long  axis  of  the  plate, 
only  half  the  plate  was  modelled  with  the  nodes  along  the  plate's 
axis  of  symmetry  being  fixed  against  asymmetrical  motions,  namely, 
translation  in  the  Y direction  and  rotation  about  the  X axis.  These 
assumptions  yield  a finite  element  model  having  957  degrees  of  freedom. 
Note  that  a refined  mesh  was  used  in  the  area  about  the  impact  point 
(shown  by  arrow).  The  plate  response  was  essentially  the  same 
whether  or  not  the  refined  mesh  was  used  (thus  demonstrating  con- 
vergence) but  better  contour  plots  resulted  when  the  finer  mesh  was 
used . 


The  transient  analysis  module  of  NASTRAN,  Rigid  Format  9,  was 
used  to  carry  out  the  analysis.  When  using  the  transient  analysis 
module  one  can  elect  to  use  either  a modal  superposition  technique 
or  a direct  integration  of  the  nodal  displacements.  The  latter 
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technique  of  direct  integration  was  adopted  for  this  study  because 
it  would  have  taken  a prohibitive  number  of  normal  modes  and  fre- 
quencies to  effectively  model  the  plate  response  for  an  impact  load 
that  had  a duration  of  55us.  A good  rule  of  thumb  when  using  the 
modal  superposition  method  is  that  both  the  integration  step  size 
and  the  period  of  the  highest  normal  mode  should  be,  at  most,  one- 
tenth  the  size  of  the  force  duration.  For  the  present  study,  this 
would  have  required  the  highest  mode  of  vibration  to  have  a period 
of  about  5ps,  i.e.,  a natural  frequency  of  2x10^  Hz.  This  fact, 
coupled  with  the  additional  fact  that  the  direct  integration  tech- 
nique is  inherently  more  accurate  since  it  handles  all  degrees  of 
freedom,  lead  to  its  choice  as  the  solution  technique. 

The  transient  analysis  module  of  NASTRAN  accepts  the  forcing 
function  in  tabular  form  where  the  load  amplitude  and  direction  versus 
time  are  input  for  each  node  in  the  area  of  the  load.  For  the  present 
study,  a single  load  in  the  normal  (Z)  direction  was  input  as  a 
function  of  time  at  the  node  designated  by  the  arrow  in  Figure  5.1 
(node  number  171).  The  load  profile  was  assumed  to  be  the  shape  of 
a half  sine  wave  (see  Figures  4.6  and  4.7)  having  a duration  of  55ys. 
By  measuring  the  height  of  the  ballistic  pendulum  at  its  release  point 
and  at  its  maximum  rebound  position  from  the  plate,  the  amplitude  of 
the  half  sine  load  can  be  obtained  by  equating  the  maximum  kinetic 
energy  to  the  maximum  potential  energy  of  the  pendulum.  For  a 
pendulum  of  mass,  m,  being  released  at  height,  hj , and  rebounding  to 
a height  of  Y\2,  equating  the  maximum  potential  and  kinetic  energies 
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yields  a change  in  velocity  of: 


Vl  - V2  = Jli  ( /h[  - /hz) 
where  g is  the  acceleration  of  gravity. 


The  impulse  of  the  force,  F(t),  acting 
for  a time,  T,  is 


I = 


/p 


T 

F (t )dt 


= m (Vj  - V2) . 


(5.1) 


(5.2) 


For  a force  having  a half  sine  wave  profile, 

F (t ) = A Sin  (-It)  (5.3) 

where  A is  the  force  amplitude. 

Utilizing  equations  (5.1)  and  (5.3)  in  equation  (5.2)  yields 
T 

I = Jk  Sin  ^ dt  = m / 2g  ( /~h[  - Jh/)  (5.4) 

0 


Carrying  out  the  integration  and  solving  for 
the  amplitude  gives: 

A TTm  / 2g  ( /hT  - ■!  hz  ) (5.5) 

2T 

-3 

The  mass  of  the  steel  ball  was  1.335x10  slugs.  From 
Figure  3.1,  it  is  seen  that  hj  = 9.75".  The  rebound  1, eight  was 
measured  by  taking  a photograph  of  the  ball's  trajectory  while 
illuminating  it  with  a high  frequency  strobe  light.  The  height,  h2 , 
was  found  to  be  .263".  Placing  these  values  into  equation  (5.5) 
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yields  an  amplitude  of  A = 230  lbf.  Hence,  the  force  on  the 
plate  is 

F(t)  = 230  Sin  (-JL £_)  (5.6) 

\55ps  / 

Equation  (5.6)  was  used  in  tabular  form  in  NASTRAN.  (This  is 
done  using  the  "DAREA"  and  "TABLED1  75"  bulk  data  cards  in  NASTRAN 
as  shown  in  the  Appendix). 

5.2  Results  of  the  Finite  Element  Analysis 

Using  the  finite  element  mesh,  boundary  conditions,  and  impact 
load  profile  described  above,  NASTRAN  computed  the  plate  displacement 
for  specified  times  after  impact.  The  output  was  in  the  form  of  dis- 
placements at  specified  nodes  plus  contour  plots  of  the  displacement 
for  the  entire  plate.  The  normal  plate  displacement  was  plotted  as  a 
function  of  the  distance  (v)  from  the  impact  point  at  X=3"  for  times 
after  impact  of  4,  6,  12,  and  18us.  For  purposes  of  comparison  with 
the  experimental  results,  the  four  plots  are  shown  in  Figures  4.10 
th  jugh  4.13  of  Section  IV  as  the  dashed  curves.  The  agreement 
between  the  experimental  and  finite  element  result  is  quite  good 
except  for  T=6ps.  This  is  felt  to  be  due  to  a stronger  than  average 
impact  load  for  the  experimental  curve.  This  agreement  is  corroborated 
by  a look  at  Figure  5.2  where  the  displacement  of  the  plate  impact 
point  is  shown  plotted  as  a function  of  time.  A computer  generated 
second  order  least  squares  curve  fit  the  finite  element  data  exactly, 
demonstrating  a parabolic  relationship  between  the  impact  point  dis- 
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Figure  5.2  - Impact  Point  Displacement  Vs  Time 
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placement  and  time.  The  experimentally  derived  points  closely 
follow  this  curve  with  the  displacement  at  T=6ps  showing  the 
largest  deviation, as  one  would  expect. 

Figures  5.3  through  5.5  show  normal  displacement  contours  of 
the  plate  response  at  times  after  impact  of  6,  8,  12,  18,  24,  and 
30ps.  The  values  of  the  contours  are  given  in  Table  5.1.  The 
trends  demonstrated  by  the  contour  plots  are  in  good  agreement 
with  the  experimental  results.  Note  in  Figure  5.5  that  reflection 
has  started  to  take  place  off  the  free  edge  for  T=30ps.  While 
some  of  the  symbols  in  the  contour  plots  may  be  difficult  to 
discern,  the  reader  can  be  aided  by  the  fact  that  contour 
numbers  1 through  29  represent  increasing  positive  displacement 
and  numbers  30  through  50  represent  increasing  negative  displace- 
ment. This  means  that  a loci  of  zero  displacement  lies  between 
contours  1 and  31,  exclusive. 
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TABLE  5.1 

j 

Normal  Displacement  Magnitude  for  NASTRAN  Contour  Plots  (See 
Figures  5.3  - 5.5) 


SYMBOL  DISPLACEMENT  (IN.)  SYMBOL  DISPLACEMENT  (IN.) 


1 

1.00E-05 

26 

9.00E-04 

> 

2 

2.00E-05 

27 

9.50E-04 

✓ 1 

3 

3.00E-05 

28 

1.00E-03 

4 

4.00E-05 

29 

1.05E-03 

5 

5.00E-05 

30 

-1.00E-05 

6 

6 . 00E-05 

31 

-2.00E-05 

7 

7.00E-05 

32 

-3.00E-05 

8 

8.00E-05 

33 

-4.00E-05 

9 

9.00E-05 

34 

-5.00E-05 

10 

1.00E-04 

35 

-6.00E-05 

11 

1.50E-04 

36 

-7.00E-05 

12 

2.00E-04 

37 

-8.00E-05 

13 

2.50E-04 

38 

-9 . 00E-05 

14 

3.00E-04 

39 

-1.00E-34 

15 

3 . 50E-04 

40 
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Figure  5.5  - NASTRAN  Contour  Plot  of  Normal  Displacement  at  T=24PS  and  T-30us 
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SECTION  VI 

DISCUSSION  AND  CONCLUSION 

The  work  described  in  the  preceding  sections  has  demonstrated 
that  double  pulsed  holographic  interferometry  can  be  effectively  used 
to  study  the  dynamic  structural  response  of  an  elastic  body  subjected 
to  impact  loading.  Within  the  limits  of  the  displacement  range  covered 
by  pulsed  laser  interferometry  using  a ruby  laser  (5-1000uin),  quanti- 
tative displacement  information  describing  a plate's  initial  response 
to  an  impact  load  was  obtained.  There  was  good  agreement  between  the 
experimental  tests  and  analytical  results  obtained  using  the  NASTRAN 
finite  element  computer  program  (Rigid  Format  9).  Some  scatter 
between  individual  tests  did  occur  and  was  felt  to  be  due  to  failure  to 
accurately  reproduce  the  impact  load  on  the  plate.  Future  work 
in  this  area  should  strive  for  a more  repeatable  impact  load.  In 
addition,  experimental  timing  measurements  should  be  increased  to  a 
sensitivity  of  +0.1  us  instead  of  the  ±0.5  us  used  in  the  present 
experimental  tests. 

The  good  experimental/analytical  agreement  provides  increased 
confidence  in  NASTRAN' s transient  analysis  module.  The  good  agreement 
does  not  come  free,  however.  For  a typical  computer  run  (see  Appen- 
dix), central  processor  (CP)  time  was  923  seconds  and  input/output 
(10)  time  was  1360  seconds  on  the  Wright-Patterson  AFB  CDC  6600 
Computer . 


Future  work  in  the  vein  of  the  current  research  will  use 
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the  existing  system  to  study  turbine  fan  blade  response  to 
impact  load.  Some  thought  is  also  being  given  to  utilizing  the 
measurement  of  the  plate  flexural  wave  speed  to  determine  elastic 
material  properties  of  composites  (say)  and,  also,  to  study  material 
damping  characteristics.  Finally,  plate  displacements,  at  times 
greater  than  those  dealt  with  in  this  report  (2-30ps),  could  be 
studied  using  the  double-pulse  capability  of  AFAPL's  ruby  laser. 
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